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Abstract. Rather weak fossil magnetic fields in the radiative core can produce the 
solar tachocline if the field is poloidal and almost horizontal in the tachocline region, 
i.e. if the field is confined within the core. This particular field geometry is shown 
to result from a shallow ( < 1 Mm) penetration of the meridional flow existing in the 
convection zone into the radiative core. I.e., two conditions are crucial for a magnetic 
tachocline theory: (i) the presence of meridional flow of a few meters per second at 
the base of the convection zone, and (ii) a magnetic diffusivity inside the tachocline 
smaller than 10 8 cm 2 s _1 . Numerical solutions for both confined poloidal fields and the 
resulting tachocline structures are presented. We find that the tachocline thickness 
runs as B p 1 with the poloidal field amplitude falling below 5% of the solar radius 
for B p > 5 mG. The resulting toroidal field amplitude inside the tachocline of about 
100 G does not depend on the B p . The hydromagnetic stability of the tachocline is only 
briefly discussed. For the hydrodynamic stability of latitudinal differential rotation we 
found that the critical 29% of the 2D theory of Watson (1981) are reduced to only 21% 
in 3D for marginal modes of about 6 Mm radial scale. 

PACS numbers: 52.30.Cv, 96.60.Hv, 97.10.Kc, 96.60.Bn, 96.60.Jw 
1. Introduction 

The tachocline is a thin shell inside the Sun where the rotation pattern changes strongly. 
Beneath the tachocline, the solar rotation is rather uniform. Above the tachocline, the 
rotation rate varies with latitude to decrease from equator to pole by about 30%. The 
transition from differential to rigid rotation detected by helioseismology (Wilson et al 
1997; Kosovichev et al 1997; Schou et al 1998) is shown in figure [TJ The tachocline 
thickness is about 5% of the solar radius, its midpoint radius is (0.692 ± O.OO5)i? , and 
it is slightly prolate in shape (Kosovichev 1996; Antia et al 1998; Charbonneau et al 
1999b). The tachocline is located mostly if not totally beneath the base of convection 
zone at R m = O.713-R (Christensen-Dalsgaard et al 1991; Basu & Antia 1997). 

The given picture of the internal solar rotation suggests that some coupling must 
exist between the base of convection zone and the inner core. Otherwise, the rotational 
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Figure 1. The internal solar rotation from helioseismology. The lines are marked by 
correspondent latitude values (courtesy NSF's National Solar Observatory). 



braking of the Sun on evolutionary time-scales would produce a rapidly rotating core 
(Dicke 1970). More than that, an even stronger link between low and high latitudes 
operates immediately below the base of convection zone in order to produce the rotation 
profile within the tachocline. The magnetic theory which explains the tachocline 
structure as a consequence of a weak internal fossil magnetic field simultaneously 
provides both the links (Riidiger & Kitchatinov 1997; MacGregor & Charbonneau 1999). 
The present paper summarizes the basic ideas of the magnetic theory and describes 
new results of numerical MHD to model the tachocline and probing the internal field 
structure. Other approaches has been recently reviewed by Gilman (2005) and Garaud 
(2007). 

Section [3] formulates the equations of the tachocline model and shows their solution 
for prescribed geometry of the fossil magnetic field within the radiative core. The role 
of meridional flow penetrating the uppermost region of the radiative core from the 
convection zone is discussed in section HI In this section, the eigenmodes of internal 
poloidal field are computed with inclusion of the penetrating flow. Tachocline models 
with consistently defined internal field modes are discussed in section [5j In the final 
section [6] a first step is presented to understand the stability problem of the solar 
tachocline. In a hydrodynamical approach the 2D Watson theory for the stability of 
latitudinal differential rotation is reformulated in 3D with surprising results. The modes 
with small radial wavelengts are most unstable so that it is indeed essential for tachocline 
theories to take into account all the radial stratifications. 

2. An heuristic approach 

The basic parameters of the tachocline theory can be estimated as follows. The estimates 
are made with a simplified shear flow model in Cartesian geometry sketched in figure [2j 
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Figure 2. A shear flow is prescribed at z = 0. Its penetration into the region of z > 
where the magnetic field Bq exists mimics the magnetic tachocline formation. 



The plane of z = mimics the bottom of the convection zone where a shear flow 
u = (U(y, z), 0, 0) imitating the differential rotation is prescribed as 

U = U ean(ky) at 2 = 0. (1) 

The Cartesian coordinates x, y and z correspond to azimuth, latitude and depth beneath 
the convection zone. An uniform (poloidal) field Bq along the y-axis is prescribed in the 
radiative zone. The shear flow penetrates the z > region and produces an x-component 
B(y, z) of the field so that 

B = (B(y,z),B Q ,0). (2) 

can be written. The equations for the stationary (toroidal) x-components of both 
magnetic and velocity fields then read 

y dy 2 dz 2 ) \i$p dy \ dy 2 dz 2 J dy 

where u is the viscosity and r\ the magnetic diffusivity. The boundary condition for 
the flow is ([1]). For the toroidal field we impose the vacuum condition B = at z = 
motivated by the very large turbulent magnetic diffusivity inside the convection zone 
compared with microscopic diffusivity of the radiative interior. The remaining two 
conditions require U and B to vanish for z — > oo. 
The solution of the equations simply reads 

U (y, z) = Uq exp (— \xkz) cos (A 2 A;z) sin ky, 

B(y, z) = y/fioP f/ Pm 1//2 exp (— Xikz) sin (X 2 kz) cos ky, (4) 
where Pm = u/rj is the magnetic Prandtl number and the parameters 
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which is the basic parameter of the whole theory. 

With Ha = 0, the equation (jSJ) gives Ai = 1, A 2 = 0, and (TjJ converts into the 
nonmagnetic solution, B = 0, U = Uoe~ kz sin ky. The resulting shear flow penetrates 
deep inside the z > region. No slender tachocline can be formed without magnetic 
fields. But with solar parameters one finds Ha ~ 10 7 Bq (Bo in G) just below the 
convection zone. Very large Ha can thus be expected. For this case, Ai = A 2 = yTfa/2, 
and the solution (TjJ provides two important consequences: 

(i) The tachocline thickness, Dtach, is strongly reduced compared to its nonmagnetic 
value, D , 

Aach = D 0^j^ ~ (V0PVV) 1/A ■ (7) 

The tachocline is so thin because its extension in ^-direction is only due to viscous stress 
while the smoothing in y-direction is performed by the much stronger Maxwell stress. 

(ii) The amplitude of the toroidal magnetic field in the steady tachocline does not 
depend on the poloidal field strength, i.e. 



B = Uo^/fioP Pm. (8) 

The ratio of magnetic to kinetic energy in the tachocline equals Pm. In terms of the 
Alfven velocity, Va = B/y/fiop, one finds 



V A = VP^ U (9) 

predicting B ~ 1000 G for the Sun. 

After ([7]) even a weak poloidal field of only B Q ~ 10~ 3 G can reduce D tach below 
5% of the solar radius. Equation ([7]) does also apply when B is not an internal field 
of the radiative core. In case the poloidal field of the solar cycle diffuses into the core 
(Forgacs-Dajka & Petrovay 2002) v and rj must be replaced by the turbulent diffusivity 
values. As their magnetic Prandtl number is not much smaller than unity the resulting 
toroidal magnetic field after ([9]) becomes much stronger than 1000 G. 



3. Tachocline model in spherical geometry 

The tachocline equations will be formulated for axial symmetry. Then the magnetic 
field can be expressed in terms of A and B in accord to 

B = e^ + Vx fe ^— ) , (10) 

where r, 9 and are spherical coordinates and e^ is the azimuthal unit vector. In the 
next section we show that meridional flow is only significant for the structure of the 
poloidal field. It can be neglected in the tachocline equations below. Here the only 
mean flow is the differential rotational with Q = Q(r, 9). 

If the poloidal field is given, the tachocline is described by the equations for the 
toroidal field B and the angular velocity fl 

r] d ( 1 d(Bsm9)\ d ( d (Br)\ _ dtt dA dVLdA 
r 89 \sin6 1 89 J 8r \ 8r J 89 8r 8r 89 ' 
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Figure 3. Radial profiles of magnetic diffusivity (fH| and viscosity (|15p calculated 
with the solar structure model by Stix & Skaley (1990). 



pis 8 ( . 3 a dVt\ 1 d ( 4 <9lV 



shr 6>— + — — r 4 p^- 



sin 9 89 \ 89 J r 2 8r \ 8r J 

_ 1 / 8A 8(Bsm9) _ . 8A 8(Br) \ 

~ ii r 2 sm 3 e \"dr 89 Sm ~89 8r J ' ^ ^ 

The rotation profile at the top boundary is produced by the convection zone for which 
the expression 

tt = 2.9 (l - 0.15 cos 2 6) /irad s" 1 at r = R in , (12) 

is used derived by Charbonneau et al (1999a). The remaining boundary conditions are 
the vacuum condition for toroidal field on the top and the conditions for the regularity 
of the fields for r = 0: 

8Q 



B\ r =R in ~ B\ r=0 - 00 



= 0. (13) 

r=0 

The used diffusivity profiles are shown in figure [3j We do not expect turbulence in 
the radiative core and apply, therefore, the microscopic magnetic diffusivity 

T) = 10 13 T~ 3/2 cmV 1 , (14) 

and the viscosity 

^ ^micro ^rad; 

z/ micro = 1.2 • 10" 16 cmV 1 , i/ rad = 2.5 • 1(T 25 — cmV 1 (15) 

P *P 

(cf. Kippenhahn & Weigert 1994) including molecular (z/ micro ) and radiative (t'rad) parts; 

k is opacity. 

The equations (fTTl) with their boundary conditions (fT2l) . (1X5]) provide tachocline 
solutions if the poloidal field is known. For a very first view the field model 

A = B p - \1 - —J sin 2 9, q>l (16) 
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is used to probe whether slender tachoclines can indeed be found. B p is the free 
amplitude of the poloidal field. Figure H] compares the rotation pattern in a radiative 




Figure 4. Isolines of the angular velocity of rotation inside the radiative core for the 
poloidal field model of (fTB)) with q = 5. The outer circle is the base of the convection 
zone. Left : B p — 0. No tachocline exists. Right: B p = 0.01 G. Tachocline exists. 

zone computed with either zero magnetic field and with Bp = 0.01 G. Bp is the maximum 
strength of poloidal field inside the core. The maximum is attained in the core center. 
The field in the tachocline region is much weaker. Even this rather weak field suffices 
to produce a very strict and slender tachocline. 

4. Magnetic field confinement by meridional flow 

Our tachocline model is not yet consistent. It works with a prescribed poloidal field. 
We have shown that the model is not very sensitive to the poloidal field amplitude. 
Even a weak field can produce tachoclinic structures. The model is, however, very 
sensitive to the field geometry. The tachocline can be found only if the field lines 
are almost horizontal near the top of radiative zone. Steady rotation and magnetic 
fields can be realized in highly conducting fluids only with constant angular velocity 
along the field lines (Ferraro 1937). Therefore, the open field geometry with field lines 
crossing the boundary between convection and radiative zones is not appropriate for 
tachocline formation. The poloidal field prescribed by (fl6|) is of a closed type. The 
field can, however, change to an open structure due to magnetic diffusion (Brun & Zahn 
2006). We shall show in the following that already a small penetration of meridional flow 
from the convection zone into the radiative core produces the confined field geometry 
required for the tachocline formation (see Kitchatinov & Rudiger 2006). If the electric 
conductivity in the core is high enough compared to that of the convection zone then 
the flow has massive consequences and the poloidal field lines in the transition zone 
become parallel to the meridional flow (cf. Mestel 1999). 
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Figure 5. Left: Penetration depth as function of viscosity of the radiative core. Right: 
Meridional velocity in the penetration region at 45° latitude for v = 1.1 • 10 9 cm 2 s _1 . 
The vertical dotted line marks the bottom of the convection zone. 



4-1. Penetration of meridional circulation 

A global poleward flow is observed at the solar surface (Komm et al 1993) persisting to a 
depth of at least 12 Mm (Zhao & Kosovichev 2004). There must be a return flow towards 
the equator somewhere deeper. Theoretical models indeed predict an equatorward flow 
of ~ 10 ms" 1 at the base of convection zone (Kitchatinov & Riidiger 1999; Miesch et 
al 2000; Rempel 2005). This flow can penetrate beneath the bottom of the convection 
zone into the radiative core (figure [5]). This penetration has been discussed recently 
in relation to dynamo models for solar activity (Nandy & Choudhuri 2002; Gilman & 
Miesch 2004; Riidiger et al 2005). 

The penetration results from the viscous drag imposed by the meridional flow at 
the base of convection zone on the fluid beneath, and this is opposed by Coriolis force. 
The thickness of such an Ekman layer results as D pen ~ ^Jv/VL (see Gilman & Miesch 

2004). The plot in the left of figure [5] is indeed approximated by D pcn ~ 2.3 ■ J u/VL or, 
equivalently, by 

D pen = 1.4 • yfO • 10 3 cm (17) 

for solar parameters^. The penetration distance is about 100 m for microscopic viscosity 
and remains shorter than 1 Mm for any reasonable value of eddy viscosity. The distance 
is small even compared to the tachocline depth, i.e. D pen -C -D tach . 

This shallow penetration is, nevertheless, highly important for the geometry of 
the internal poloidal field. The ratio of the time r diS = D^ cn /r] for the magnetic 
field diffusion across the penetration layer to the characteristic time of induction of a 
meridional magnetic field from a radial one, r s = D pcn /u m , gives the magnetic Reynolds 

t -Dpcn of figure [5] was defined as distance from the convection zone bottom to the position where the 
meridional flow falls to zero 
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Rm=^ = ^^^.10 6 (18) 



number 



(with u m ~ 10 ms _1 for the meridional velocity). This Reynolds number is large, 
Rm~ 10 3 , for microscopic diffusivities but it remains above this limit with eddy 
diffusivities up to 10 6 cm 2 s _1 . The large ratio of fflBl means that the latitudinal field 
inside the penetration layer is large compared to the radial field component. With other 
words, the field has just the confined geometry required for the tachocline formation. 
Only with eddy diffusivities of 10 12 cm 2 s _1 the Reynolds number ({TBI sinks to unity, so 
that only in this case the influence of penetration on the internal field geometry is too 
weak. 

The ratio of r^g to the advection time r a d v = Rin/u m 



Tdiflf -Dpcn 



Rm ~ 0.03 Pm < 1, (19) 



^adv R'ir 

is small independently of whether microscopic or eddy diffusivity is used. If the 
tachocline is stable in the hydrodynamic regime (see section 6) the microscopic 
diffusivities should be used. The small ratio ([19]) justifies the neglect of the meridional 
flow in the tachocline equations (fTTj) . The diffusion time, r difT , is too short as could the 
penetration layer be dynamo- relevant. The configuration of the internal poloidal field is 
probably the only process for which the penetration of the meridional circulation into 
the stable radiative zone is important. 

4-2. A model of the poloidal field 

For both axisymmetric flows and fields the poloidal magnetic field equation decouples 
from the tachocline equations (TTTj) and reads 

9A__^dA_ u 2A + T) ^A + IL s i n e—(——\ (20) 
8t r 86 r dr dr 2 r 2 86 \ sin 6 86 J ' 

where u is the velocity field. The vacuum boundary condition for the poloidal field can 
be applied on the top of the radiative zone because of the very large turbulent diffusivity 
in the convection zone, 

8A (8A\ , . 

\ / vac 

The condition is usually formulated in terms of a Legendre polynomials expansion 

oo 

A{r ,6) = sin0£4»(r)P,i(cos0), 




n=l 

sin# 00 



r n=l 



J2nA n (r)P^ (cos 6). (22) 



The other condition is A = at r = 0. 
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The velocity u can be written with its stream function tp as 

y pr z sin u OU pr sin 6 or J 

The meridional flow in the bulk of the radiative core is very slow. The characteristic 
time of the Eddington-Sweet circulation even exceeds the solar age (Tassoul 2000). Only 
the flow penetrating from convection zone is significant for the tachocline. The stream 
function, xjj(r,9), of the penetrating flow drops to zero at a small depth of order D pen 
inside the core. The function is, however, finite at the top boundary of the core where 
it scales as 

i{;(R in ,e)=u m pR- m D pe J(e), (24) 

where ip is dimensionless function of order unity. u m is the meridional velocity amplitude 
at the boundary. 

The penetration depth is so small even compared to tachocline thickness that we 
are motivated to integrate equation ( 120]) across the penetration layer instead of resolving 
this layer explicitely. Such an integration reformulates the top boundary condition which 
now reads 

S£ -*•(£) • 



where Rm is the magnetic Reynolds number ( 118]) and the RHS is defined in ( 122]) . For 
Rm = (123]) reduces to the vacuum condition (12T]) . The true Reynolds number is 
Rm ~ 10 3 . 

Penetrating meridional flow is now included via the condition (]25p and poloidal 
field structure can be found by solving eigenvalue problem for diffusion equation, 

A d 2 A v n d ( 1 dA\ . , 

-T =v w + 7* smd w{^rew)> (26) 

where eigenvalue r is the time of resistive decay of normal modes. The results discussed 
below were obtained in computations with the simplest stream function, 

ip(9) = -sin 9 P£ (cos 9) , (27) 

where P\ is the normalized Legendre polynomial. 

The penetrating flow is expected to confine the internal poloidal field within the 
radiative core. The parameter 

(28) 



max 


A (r, 9) 


r=R b 


max 


A (r, 9) 


r<R b 



is used to estimate the field confinement. The 50-parameter ( 128]) estimates the ratio of 
the magnetic flux through the surface of the core to the characteristic value of the flux 
within the core. 

The linear equations (126]) and f]25]) define the field structure but not its amplitude. 
The amplitude B p of the poloidal field is a free parameter of the model. 
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4-3. Normal modes of the internal field 

How the structure of the internal field changes towards a confined geometry as Rm 
increases is shown in figure EJ The internal field computed with Rm = has an open 
structure. Even a moderate flow with Rm = 10 changes the field considerably towards 
the confined geometry. The internal field for the solar value Rm = 1000 shown in 
figure [6] is almost totally confined. Less than 1% of magnetic flux belongs to 'open' 
field lines in this case (figure EJ) . The figures [6] and [7] only display the most slowly 
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Figure 6. Field lines for the longest-living dipolar modes of poloidal field for different 
Rm. The Reynolds number varies as Rm = 0, 10, fOOO from left to right. The outer 
circle is the base of convection zone. The decay times in Gyr are marked at the tops. 
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Figure 7. Confinement parameter l|28p for the longest-living dipolar modes as function 
of the magnetic Reynolds number Rm. 



decaying dipolar modes. Other normal modes have shorter scales in radius or latitude 
and also shorter decay times. The three modes following the longest-living dipole ordered 
for decreasing lifetimes are shown in figure El The higher-order modes also have a 
confined structure. The decay times of the given normal modes are long enough to 
allow tachocline computations with steady poloidal fields. 

For sufficiently high Rm the direction of the penetrating meridional flow does not 
play a role. The flow from equator to poles also makes efficient confinements of the 
internal field. The effect can be understood as a magnetic field expulsion from the region 
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Figure 8. Normal modes of internal field following the dipolar mode of figure [5] for 
shorter lifetimes. The decay times (in Gyr) are shown on the top. The confinement 
parameter (gHJ) varies as 0.0051, 0.0073, 0.0026 from left to right. Rm = 1000. 
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Figure 9. Angular velocity isolines in the radiative core computed with dipolar 
internal fields for Rm = 1000 for various poloidal field amplitudes (marked at the 
top). 



of circulating motion (Weiss 1966). The transition from convection zone to radiative core 
was treated as a sharp boundary in our model. A smooth decrease of eddy diffusivity 
in convection zone towards its base may also contribute to the field confinement via the 
diamagnetic effect of inhomogeneous turbulence (cf. Krause & Radler 1980). 

5. Tachocline models 

Figure M shows the angular velocity distributions inside the core computed with dipolar 
eigenmodes of internal field of various strengths. Dependence of the tachocline thickness 
on poloidal field amplitude is shown in figure [10] (-Dtach was defined as depth of 
exponential decrease of equator-to-pole difference in angular velocity). Even a very 
weak field, B p ~ 10~ 3 G, can produce tachocline. The depth -D tac h drops below 5% of 
the solar radius for B p > 5 mG. 

The tachocline thickness in relation to the amplitude to the prescribed fossil field 
is given in figure [Hfl It is very close to the rough estimate (171). The simulations also 
confirm the above finding that the toroidal field amplitude, B t = 100. ..200 G, hardly 
varies while the poloidal field amplitude changes by several orders of magnitude. 
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Figure 10. The fractional tachocline thickness as function of the amplitude, B p , of 
the internal poloidal field. Rm= 1000. 



Also the higher order modes given in figure [8] also produce tachoclinic structures 
provided that the fields have the confined structure which is always the case for Rm= 100 
or larger. 

Figures [9] and [10] are valid for Rm = 1000. Already Rm = 100 suffices for tachocline 
formation, but for smaller Rm the poloidal field remains 'too open' for the formation of 
a tachocline. With other words, a meridional flow with minimum 1 ms _1 amplitude at 
the base of convection zone is necessary for the success of the magnetic tachocline theory. 
The rather shallow penetration of the meridional flow of the convection zone into the 
radiative core influences the internal field geometry strongly enough and in such a way 
that the field becomes appropriate for the tachocline formation. We suggest that a layer 
thiner than 1 Mm beneath the convection zone where a meridional flow of 1-10 ms _1 
enters the radiative core is responsible for the existence of the solar tachocline. 

6. Hydrodynamic stability: the Watson approach in 3D 

Latitudinal differential rotation can be unstable even without magnetic field if the 
shear dQ/dd is positive and sufficiently strong (Watson 1981). The critical value of 
29% latitudinal shear found by Watson for the nonaxisymmetric mode with m = 1 
resulted from a theory with strong stratification but without any radial velocity and 
radial wavelengths. The value has also appeared in a 3D numerical examination of 
marginal stability of a shell rotating fast enough with the rotation law 



but of fully incompressible material (Arlt et al 2007). The critical shear increases to 
higher values, however, if the real rotation law (including its radial variations) of the 
solar tachocline is adopted. In this case the superrotation observed in the equatorial 
region of the rotation law strongly stabilizes the shear instability. 

In this section the Watson approach with strong density stratification (for m = 1) 
is extended to 3D, i.e. to the inclusion of finite radial velocities and finite radial 




(29) 
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Figure 11. The buoyancy frequency (150)) in the upper radiative core after the solar 
model of Stix & Skaley (1990). 



wavelengths. The stabilizing effect of a subadiabatic stratification on the mean flow 
is characterized by the buoyancy frequency, 
9 OS 



N = — 

C p dr 



(30) 



where S = C v In (P/p 1 ) is the specific entropy of ideal gas. The frequency is very large in 
the solar radiative core (see figure [TTT) . The larger the N the more the radial fluctuations 
are suppressed by the buoyancy force. Radial velocities should therefore be small. Our 
stability analysis is local in the radial dimension, i.e. we use Fourier modes exp ikr with 
fcr>l. The equation system remains, however, global in horizontal dimensions. 

Instabilities of rotating shear flows are fast. Their growth rates are of order Q 
or smaller. In the short-wave approximation the velocity field can be assumed as 
divergence-free, divu' = 0. The next assumption concerns the pressure. Local thermal 
disturbances occur at constant pressure so that p'/p = —T'/T or S' = —C p p'/p where 
perturbations have been marked by dashes. This assumption is again justified by the 
incompressible nature of the fluid. We start from the linearized equations of velocity, 
i.e. 



and entropy 



du' 
~dt 



dS' 



+ (u • V) u' + («' • V) u 




+ uAu' 



+ u ■ VS' + u ■ V S 



AT'. 



(31) 



(32) 



dt T 
The basic flow is the rotation law ( 1291) . 

Perturbations of velocity and pressure are now expressed in terms of scalar 
potentials like 

1 dT d 2 P \ e s ( dT 1 d 2 P 

~d9 



u 



-LP 




sin 9 drdcf) ) 



(33) 
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with the operator 



sm6d9 Sm 89 + sin 2 9 deb 2 ' { } 



The identities 

r (r • V x u') = LT, r 3 (r ■ V x V x u') = [L + r 2 |- ) LP, (35) 



are used to reformulate the relations in terms of the potentials. 

The perturbations are considered as Fourier modes for time, azimuth and radius 
in the form exp (i(— out + mcf) + kr)). For an instability the eigenvalue uo must possess a 
positive imaginary part. Only the highest order terms in kr for the same variable have 
been considered. The pressure term in ( 151]) can be formulated as 

r ■ V x V x ( -Vp ) = - r ■ V x ( -^(Vp) x (Vp) ] = 



JL . V x f i(V5) x (Vp)) = - ■ V x (g x VS") = -^Is". (36) 

In order to denormalize the variables time is measured in units of 1 and velocities 
are scaled with rflo. The remaining variables are 

v -^ p ' W = ^ T < ' = fl = £ (37 » 

Then the equation for the poloidal flow reads 

= " if (^) " 

-2^(W)-2(1-^)^M^-2A + 

V / \ 7 Ojl Ojl OfJ, 

+ mh (LV) + 2m^-V + 2m (l - p 2 ) (38) 
v ' dp v ' op op 

with /i = cos# and the normalized radial wavelength 

The first term on the RHS describes the stabilizing effect of the stratification. It vanishes 
for small A. Apart from this stabilizing buoyancy term, the wavelength only exists in 
diffusive terms. The second term of the RHS includes the action of viscosity. Here we 
have used the quantities 

vN 2 X N 2 
€u ^r 2 ' 6x Q 3 r 2 1 } 

for the diffusive parameters v and x- Now the limit N — > can no longer be used for 
the transition to nonstratified but viscous media. The numerical values for the solar 
tachocline are 

e v = 2 ■ 1(T 10 , e x = 10 -4 . (41) 
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The second and the following lines in (I3"8j) describe the influences of the basic rotation. 
Note that only latitudinal derivatives of Q appear. 

The complete system of equations also includes the equation for toroidal flow, 

d 2 

(lw) = -i^ (LW) + m(l (LW) - mlf — ((l - ^i 2 ) Ct) + 





and the entropy equation 

us = -\^-s + mVts + LV. (43) 
A 2 

Note that for the Reynolds number Re = fir 2 jv follows 
N 2 /fl 2 

Re = ilZ!i (44) 
£u 

which with (I4ip is very large. 

For positive a in (I29p the modes Al and S2 become unstable for sufficiently large 
a = fl crit- Figure [121 shows the dependence of the a crit on the normalized wavelength A. 
Of course, for large enough radial wavelength the 29%-value of the Watson theory with 
2D fluctuations is reproduced. It is reduced, however, to a = 0.21 in the full 3D model 
calculations. We see that in contrast to the original Watson approach short rather than 
long radial scales are preferred. The minimum a cr it appears for A ~ 0.6, so that the 
characteristic wavelength of A ~ 6 Mm results for the solar tachocline. For shorter 
scales the instability for m = 1 disappears due to the action of the dissipation processes 
but modes with higher m remain still unstable. 



7. Concluding remarks 

The presence of a meridional flow with an amplitude of a few meters per second at 
the bottom of the convection zone is shown as necessary for the magnetic tachocline 
model. The flow provides the confined geometry (with field lines parallel to the outer 
spherical boundary) of the internal magnetic field which is necessary for the tachocline 
formation. The bottom flow is also a key ingredient of advection-dominated dynamo 
models for the solar cycle (see Riidiger & Hollerbach 2004 for detailed references) . The 
flow was predicted theoretically but its existence is not yet confirmed by observations. 
We hope, of course, that helioseismology will probe the deep meridional flow soon. 

The polar cusps in the tachoclines of figure [9] are consequences of the assumed 
axisymmetry. This assumption may not be realistic. The penetrating meridional flow 
is, however, expected to confine even nonaxisymmetric fields to the core so that the 
fields should also be able to produce the tachoclinic structures. 

The tachocline computations result in toroidal fields of 100 to 200 G, much stronger 
than the original poloidal fields. These toroidal fields can be subject to current- 
driven pinch-type instabilities (Tayler 1973). Hence, the question is whether the 
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Figure 12. Stability map for latitudinal differential rotation (3D hydrodynamics). 
Most unstable are the perturbation modes with the vertical scale A ~ 0.6. The critical 
magnitude of latitudinal shear is reduced to 0.21 compared to the 0.29 value for large 
wavelengths. For short wavelengths the instability for m = 1 rapidly disappears by 
dissipation effects. 

steady axisymmetric solutions of our model are stable against all the nonaxisymmetric 
disturbances. The threshold field strength for the instability were estimated to be of 
order 10 2 - 3 G (Spruit 1999; Arlt et al 2007; Kitchatinov & Rudiger 2007) what is (only) 
somewhat larger than the field amplitudes in our model. 

It is important for the magnetic tachocline model that the magnetic Reynolds 
number (fl8l) is large. This is true for microscopic diffusivity or for not too large eddy 
diffusivities up to about 10 8 cm 2 s _1 . The tachocline should, therefore, be stable or 
only mildly turbulent to allow poloidal field confinement by meridional flow. Some low 
level of turbulence may, however, be even necessary. The magnetic field is so efficient 
in producing tachocline structures that poloidal fields of only 1 G lead to radial scales 
smaller than 1% of the solar radius. It is, however, -Dtach — 0.04R Q . Hence, the internal 
poloidal field must be (much) smaller than 1 G (Kitchatinov et al 2001) or, if it is 
not, some kind of instability may prevent the tachocline thickness to reduce below the 
helioseismologically detected level. 

The stability issue is a clear perspective for future tachocline studies. As a 
corresponding application we have developed in the last section the 2D theory of Watson 
(1981) for the hydrodynamic instability of latitudinal differential rotation to a fully 3D 
theory. Now the critical wavelength of the unstable mode with m = 1 is only 6 Mm 
while the critical latitudinal shear (see (129]) ) is reduced from 29% to 21%. This (linear) 
theory works with a Prandtl number of order 10~ 6 and a very high Reynolds number of 
order 10 12 . 
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